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ABSTRACT 

We discuss the time variation of dust-to-gas mass ratio in spiral galaxies 
using the multi-phase model of interstellar medium. The typical timescale of 
the phase change of an interstellar gas is ~ 10 7 ~ 8 yr in spiral galaxies. Since the 
phase transition changes the filling factor of the cold gas where the dust growth 
occurs, the dust growth rate varies on that timescale. In order to examine the 
response of the dust-to-gas ratio to the phase transition, we construct a model 
of the time evolution of the dust-to-gas ratio. We adopt the three phase model 
for the interstellar gas and the Ikeuchi-Tomita model for the mass exchange 
among the phases. According to the model, three types of solutions are possible: 
[1] all the gas is transformed to a hot gas; [2] a stable stationary state of three 
phases is realized; [3] the filling factor of each phase cyclically changes. For 
each of the three types of solutions, the dust-to-gas ratio behaves as follows: 
[1] almost all the dust is destroyed (the dust-to-gas ratio becomes ~ 0); [2] the 
dust-to-gas ratio converges to a stationary state; [3] the dust-to-gas ratio varies 
cyclically in response to the phase transition. In the case of [3], the amplitude 
of the variation of the dust-to-gas ratio is large (nearly an order of magnitude) 
if the dust growth timescale is shorter than the phase transition timescale. This 
condition is easily satisfied in spiral galaxies. However, it is difficult for dwarf 
galaxies to realize the condition because their small metallicity makes the dust 
growth timescale long. 

Subject headings: galaxies: evolution — galaxies: ISM — galaxies: spiral 
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1. INTRODUCTION 

Some of the recent chemical evolution models of galaxies have successfully included 
the dust content (Wang 1991; Lisenfeld & Ferrara 1998; Dwek 1998; Hirashita 1999a, 
hereafter H99). Takagi, Arimoto, & Vansevicius (1999) modeled the dust extinction in 
ultraviolet (UV)-optical wavelengths and the reprocessed light in the far-infrared (FIR) 
region, presenting the spectral evolution in the UV-FIR wavelengths. 

As to the dust formation, supernovae (SNe) are known to be one of the most dominant 
sources (Barlow 1978; Draine & Salpeter 1979; Dwek & Scalo 1980), and SN shocks destroy 
grains (McKee et al. 1987; Jones et al. 1994; Borkowski & Dwek 1995). Thus, dust content 
is connected with star formation histories. In our previous work, H99, the dust-to-gas 
mass ratio was expressed as a non-linear function of metallicity, which is also related to 
star formation histories (see also Lisenfeld & Ferrara 1998). H99 confirmed the suggestion 
proposed by Dwek (1998) that the accretion process onto preexisting dust grains is efficient 
in spiral galaxies. This is also supported by the observations of elemental depletions in 
dense clouds (e.g, Draine 1990). Since the accretion is effective in cold clouds, the global 
efficiency of the accretion depends on the fractional mass of the cold gas (Seab 1987; McKee 
1989; Draine 1990; Hirashita 1999b). Thus, the efficiency varies on a timescale of the phase 
transition of the interstellar medium (ISM; ~ 10 7 ~ 8 yr; Ikeuchi 1988; McKee 1989). 

McKee & Ostriker (1977) constructed a standard model for the ISM composed of three 
phases: the hot rarefied gas (T ~ 10 6 K and n ~ 10~ 3 cm -3 ), the warm gas (T ~ 10 4 
K and n ~ 10 _1 cm -3 ), and the cold cloud (T ~ 10 2 K and n ~ 10 cm -3 ). One of the 
simplest models for the mass exchange among the three ISM phases is described by Ikeuchi 
& Tomita (1983, hereafter IT83; see also Habe, Ikeuchi & Tanaka 1981). They included the 
following three processes: [1] the sweeping of a warm gas into a cold component, [2] the 
evaporation of cold clouds embedded in a hot gas, and [3] the radiative cooling of a hot gas 
by collisions with a warm gas. Mass fraction of each component can show an oscillatory 
behavior, which is known as a limit-cycle (Nicolis & Prigogine 1977; see also Fujimoto & 
Ikeuchi 1984 for another modeling of interstellar gas). The limit-cycle behavior is supported 
by Kamaya & Takeuchi (1997, hereafter KT97), who interpreted the various level of the 
star formation activities in spiral galaxies shown by Tomita, Tomita, & Saito (1996) using 
the limit-cycle model. 

In the limit-cycle model by IT83, mass fraction of each phase oscillates continuously 
because of the mass exchange among the three components of the ISM (see also Korchagin, 
Ryabtsev, & Vorobyov 1994). The timescale of the phase transition in the model is 
determined by three timescales of the above physical processes [l]-[3]. Actually, a static 
solution as well as the limit-cycle solution for the mass fractions is possible for some range 
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in the parameter space (this case is also examined in this paper). However, the oscillatory 
behavior (i.e., the limit-cycle model) of the fractional masses is supported observationally. 
Indeed, the observed scatter of the FIR-to-optical flux ratios of spiral galaxies (Tomita, 
Tomita, & Saito 1996) is interpreted through the limit-cycle model in KT97, who suggested 
that the fractional mass of the cold phase changes in the range of 0.1 to 0.7 (or more) on the 
timescale of 10 7 ~ 8 yr and that this leads to the time variation of star formation activities if 
the Schmidt law (Schmidt 1959) is applicable. 

Hirashita (1999b) combined the framework of H99 with a theoretical work on 
multiphase ISM and suggested the time variation of the dust-to-gas ratio by the phase 
transition. We aim at constructing a model of the time evolution of dust-to-gas ratio 
including the effect of the ISM phase transition. The model in this paper enables us to 
simulate the time variation of the dust-to-gas ratio in a galaxy. This paper is organized as 
follows. First of all, in the next section, a model for the time variation of the dust-to-gas 
ratio is presented. Next, we show the numerical solution of the model in §3. In §4, we 
discuss the solution and present observational implications. Finally, §5 is devoted to 
summary. 



2. MODELING OF DUST-TO-GAS RATIO AND ISM PHASE 

TRANSITION 

2.1. Dust-to-Gas Ratio 

First, we prepare a model for the time variation of the dust-to-gas ratio in a galaxy. 
The model is developed by H99 (see also Lisenfeld & Ferrara 1998 and Dwek 1998) in 
order to account for the relation between metallicity and dust-to-gas ratio of nearby spiral 
galaxies. He adopted a simple one- zone model (i.e., the galaxy is treated as a closed system 
and spatial variations of physical quantities within the galaxy is neglected). The model 
equations describe the changing rate of th mass of gas, metal, and dust: 

^ = -X,i- + E„ (2) 

at r acc r SN 

Here, M g is the mass of gas; Mj and M dj j denote the total mass of the metal % {% — O, C, 
Si, Mg, Fe, etc.) in both gas and dust phases, and only in the dust phase, respectively. 
The star formation rate is denoted by ip; E is the total injection rate of mass from stars; 
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Xi is the mass fraction of the element % (i.e., Xi = Mj/M g ); i% is the total injection rate 
of element z from stars; /j is the mass fraction of the element % locked up in dust (i.e., 
fi = M^i/Mj). The meanings of the other parameters in the above equations are as follows: 
/in 5 i represents the dust mass fraction in the injected material (in other words, the dust 
condensation efficiency in the ejecta such as stellar winds or SNe); a refers to the efficiency 
of dust destruction during a star formation [a = 1 corresponds to destruction of only 
the dust incorporated into the star, and a > 1 (a < 1) corresponds to a net destruction 
(formation) in the star formation]; r acc is the accretion timescale of the element % onto 
preexisting dust grains in molecular clouds; r SN is the timescale of dust destruction by SN 
shocks. We hereafter assume that only dust incorporated into stars is destroyed when the 
stars form; i.e., a = 1. The mass loss from the galaxy is neglected in this paper, since 
we mainly treat spiral galaxies and expect that outflows are prevented by their deep dark 
matter potentials (Habe & Ikeuchi 1980). When we consider star burst phase (Habe & 
Ikeuchi 1980) or dwarf galaxies (Mac Low & Ferrara 1999), the outflows may be important 
for the dust content as well as gas and metal content (Lisenfeld & Ferrara 1998). 

We adopt the instantaneous recycling approximation (i.e., a star with the mass larger 
than mi dies instantaneously after its birth, leaving the remnant of mass w m ). Once the 
initial mass function is fixed, we can calculate the returned fraction of the mass from formed 
stars (71) and the mass fraction of the element i that is newly produced and ejected by stars 
(yi)- Using 71 and y^ E and E^ in equations ([j])-© is written as (see H99 for the details) 

E = TtiP, (4) 
E % = (TZXi + y^. (5) 

After some algebra, equations (0) and (^) becomes 

TSF^f = y, (6) 

tsf^ = f^iinXi+yd-la-l + K-P^l-tt+PssYDi, (7) 

where r SF = M g /ip (timescale of gas consumption; e.g., Roberts 1963), Pj = M dii /M g = fiX^ 
j3 acc = TsF/T acc , and /?sn = t sf/tsn- The timescales above are estimated in the case of 
typical spiral galaxies as t$f ~ HO Gyr (Kennicutt, Tamblyn, & Congdon 1994), and 
r aC c ~ r SN ~ 10 s yr (Draine 1990). 

The two timescales, r acc and tsn depend on the phase of the ISM where the dust exists 
(Draine 1990; Hirashita 1999b). The accretion timescale, r acc , is shorter in the cold phase of 
the ISM than in any other phase, since efficient accretion needs a dense environment. If we 
express the growth timescale of a dust grain in the cold component as r grow , the accretion 
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timescale is expressed as r acc = T grow /X cold , where X cold represents the mass fraction of the 
cold phase to the total mass of ISM (Hirashita 1999b). This is equivalent to 

acc ^cold/^grow; (8) 

where /3 grow = T SF /r grow . According to Draine (1990), r grow ~ 5 x 10 7 yr in the cold gas. 
The destruction timescale of the dust /?sn is kept constant for simplicity, since we would 
like to concentrate on the time variation of the dust growth efficiency to examine the 
result in Hirashita (1999b). The star formation history should be taken into account in 
order to model the time dependence of /?sn- Moreover, the time variation of the SN rate 
also contributes to the time variation of the fractional masses the three ISM components. 
Because of such complexities in the time dependence of /?sn, we simply treat /3sn as 
constant. 



2.2. ISM Phase Transition 

Next, we review the model by IT83. This model is used to calculate the time evolution 
of the filling factors of the three ISM phases (see also §1). The result is used to calculate 
the dust formation efficiency through equation (^). 

The interstellar medium is assumed to consist of three components (McKee & Ostriker 
1977); the hot rarefied gas (T ~ 10 6 K, n ~ 10 -3 cm -3 ), the warm gas (T ~ 10 4 K, 
n ~ 10 -1 cm -3 ), and the cold cloud (T ~ 10 2 K, n ~ 10 cm -3 ). The fractional masses of 
the three components are Xhot) X warm , and X co id, respectively. A trivial relation is 

Xhot + X warm + X co id = 1- (9) 

The following three processes are considered (see IT83 for the details): [1] the sweeping 
of a warm gas into a cold component at the rate of aX warm (a ~ 5 x 10 -8 yr -1 ); [2] the 
evaporation of cold clouds embedded in a hot gas at the rate of 6A co idA^ ot (b ~ 10~ 7 -10 -8 
yr -1 ); [3] the radiative cooling of a hot gas through collisions with a warm gas at the rate 
of cX warm Ahot (c ~ 10~ 6 -10~ 7 yr -1 ). Writing down the rate equations and using equation 
(|), we obtain 

BX co \ d Xl ot + A(l — X cold — X hot ), (10) 
Xhot(l — X co id — Xhot) + -BX co idX^ ot , (11) 



dX, 



cold 



dr 
dX hot 

dr 



where r = ct, A = a/c, and B = b/c. 
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The solutions of equations (|T0|) and ([□]) are classified into the following three types 
(IT83): 

[1] A > 1; all the orbits in the (X co id, Ahot)-plane reduce to the node (0, 1), 

[2] A < 1 and B > B cv ; all the orbits reduce to a stable focus [(1 - A)/(AB + 1), A], 

[3] A < 1 and B < B CT ; all the orbits converge on a limit-cycle orbit, 

where B cr = (1 — 2 A) /A 2 . As a representative parameter sets for [1] and [2], we examine 
(A, B) = (1.2, 3.0) and (A, B) = (0.4, 3.0), respectively (denoted as Models 1 and 2 
in Table 1, respectively). We especially concentrate on the case [3], since we treat spiral 
galaxies with cold gas and with various levels of star formation activities (e.g., Tomita et al. 
1996). The two properties of spiral galaxies are naturally realized (KT97) if the filling factor 
of the cold gas shows the limit-cycle behavior and the star formation rate changes according 
to the Schmidt law (tp oc X^ old ; n = 1-2; Schmidt 1959). As representative parameter sets 
of the limit-cycle orbit, we examine the case of (A, B) = (0.3, 3.0) and (A, B) = (0.3, 0.5) 
(denoted by Model 3 and 4 in Table 1, respectively) Q. The cycle periods of the former case 
is about three times shorter than that of the latter. 

In Table 1, the adopted parameters are summarized. The column "type" means the 
behavior of the solutions and named "runaway," "stationary" and "limit-cycle" for the 
cases [1], [2] and [3], respectively. We note that all the parameter sets above are examined 
in IT83 (their Figures 1 and 3). 



2.3. Summary of the Basic Equations 

Here, we summarize the equations and the parameters. 

First, the filling factor of the cold gas is calculated by the following set of equations: 

— ; — = —BX co \ d Xl ot + A(l — X cold — X hot ), (12) 

GST 

dXh t 

^ = — Ahot(l — A co id — Ahot) + -BA co idA hot , (13) 

with parameters A and B listed in Table 1. According to these equations, the time 
variation of the cold-gas fractional mass A co id is calculated to obtain the efficiency of the 
heavy-element accretion onto dust grains (eq. ||). 



1 Two cases for each model are examined for different values of c as explained in §2.3 (Model 3a, 3b, 4a, 
and 4b). 
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Next, the time variation of the dust-to-gas ratio is simultaneously solved. Since 
we are interested in a timescale much shorter than tsf, which is also interpreted as 
metal-production timescale through equation (|6|), we treat X$ as constants. Then, equation 
(0) is solved for oxygen abundance (i — O, because we adopted the oxygen as a tracer 
element in the previous work, H99: see also Lisenfeld & Ferrara 1998): 

dry 

C7 SF^ = /in, (^0 + ^o) ~[n- /Ml - fo) + MVO, (14) 

where we express the time in a nondimensional manner by using r = ct and we assume 
a = 1 (equivalent to the assumption that only dust incorporated into stars is destroyed 
in the star formation processes; §2.1). We note the relation fo = T>q/Xq. We fix the 
parameters as r SF = 3 x 10 9 yr (e.g., Kennicutt et al. 1994), / in = 0.05 (H99), 7Z = 0.79, 
3^o = 1-8 x 1CT 2 , /9g r ow = 100 (corresponding to r grow = tsf/100 = 3 x 10 7 yr; we note 
that /? g row and /3 acc are related through eq. ||), /?sn = 10 (H99), and Xo = 0.013 (the solar 
system abundance). The dependence on these parameters above are described in Lisenfeld 
& Ferrara (1998), H99, and Hirashita (1999c) except for c and rgp. The parameter c scales 
the phase transition timescale, and we examine the cases of ctsf = 100 and 1000 (Table 1: 
cT l = 3 x 10 6 -3 x 10 7 yr; IT83). We note that c appears in the form of the nondimensional 
parameter ctsf- We examine only the case of ctsf = 100 in Model 1 and 2, since the 
qualitative behavior of the dust-to-gas ratio is the same, irrespective of which value of ct$f 
is adopted. 

Finally, the total dust-to-gas ratio D is calculated from the following equation: 

V = CV , (15) 
where we assume C = 2.2 (the Galactic value; Whittet 1992, p. 52). 



We solve the basic equations Ql2)— (|15) for the parameters listed in Table 1. The initial 
condition for the filling factors are set according to IT83: (X co id, -Xhot) = (0-0, 0.1) except 
for the Model 4 and (X cold , X hot ) = (0.0, 0.7) for Model 4. The initial dust-to-gas ratio is 
set as T> = 0.0080 (Galactic value; Whittet 1992) for all the models. Models 3 and 4 are 
examined for two different values of c listed in Table 1. 



3. RESULTS 

First, the result of Model 1 is displayed in Figure 1. Because the cold gas disappears, 
the dust growth becomes ineffective. Then the dust destruction makes the dust-to-gas ratio 
converge to ~ in the dust destruction timescale (~ 3 x 10 8 yr, i.e., 10 in the unit of c _1 ). 
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This model represents the efficient sweeping of warm gas in comparison with the radiative 
cooling of the hot gas. Thus, the swept gas is first converted into the cold component and 
then into hot gas (IT83). If the hot gas escapes out of the galaxy owing to its large thermal 
energy, this model corresponds to the galactic wind model of elliptical galaxies (IT83; 
Arimoto & Yoshii 1987). 

The final value of the dust-to-gas ratio in Model 1 is calculated by considering a 
stationary state, dVo/dt = 0, in equation (|HD: 

^V 2 + {K- /3 acc + /? SN )£> - fin,o(KX + y Q ) = 0. (16) 

Because of the constant dust formation from metal injection at the rate of fm,o(R.X + yo), 
the dust-to-gas ratio is not exactly even after a long time. Since X cold = 0, /3 acc = 0. 
Thus, the solution of equation (|16|) is 

V = CV Q = C ^ nX ^ + yo) = 3 x 1(J -4 (17) 

rl + PSN 

for the adopted parameters Xo = 0.013, /; n o = 0.05, 1Z = 0.79, = 0.018, and /?sn = 10. 
This indeed matches the final value of the dust-to-gas ratio in Model 1. 

Next, the result of Model 2 is presented in Figure 2. As shown in IT83, the solution 
for X cold reduces to a stable point (1 — A)/(AB + 1) = 0.27. The stable stationary state is 
realized at a finite dust-to-gas ratio, which is determined so that the dust formation rate 
is balanced with the dust destruction rate. The value is calculated from equation flU]). 
Putting the adopted parameters, the positive solution is V = CV Q ~ 0.017, which matches 
the stationary value in Figure 2. 

Finally, the limit-cycle case is examined in Figure 3 and 4. The four sets of parameters 
listed in Table 1 are examined (Models 3a, 3b, 4a, and 4b). The resulting time evolution 
of the cold component of ISM and the dust-to-gas ratio is presented in Figures 3a, 3b, 
4a, and 4b for the four cases. As expected, the dust-to-gas ratio oscillates in response to 
the fractional mass of the cold component, though the amplitude of the variation largely 
depends on the parameters. Comparing the Figures 3a and 4a (or Figures 3b and 4b), we 
see that the amplitude becomes larger if the period of the oscillation is longer. This is 
because the dust has enough time to grow if the period is long. Also from the comparison 
between the Figures 3a and 3b (or Figures 4a and 4b), we find that the smaller c makes the 
amplitude larger. Since small c means that a timestep of the phase transition model (eqs. 
|T0| and [11]) is long in the real time unit (r = 1 corresponds to t = c _1 ), the dust has enough 
time to grow in the case of smaller c. The timescale of the dust growth relative to that of 
the phase transition determines the amplitude. This point is again discussed in the next 
section. 
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4. DISCUSSIONS 
4.1. Timescales 

In the runaway and stationary solutions (Models 1 and 2), the dust-to-gas ratio is 
settled into a stationary value determined by equation (|16D after the dust formation and 
destruction timescales (~ 10 8 yr). In the limit-cycle type solutions (Models 3 and 4), 
the amplitude of the dust-to-gas oscillation largely depends on the timescale of the phase 
transition. In this subsection, we discuss the phase-transition and dust-formation timescales 
in the limit-cycle model. 

From the solution of the model equations in the previous section (Figs. 3 and 4), 
we find that the dust-to-gas ratio oscillates in response to the filling factor of the cold 
component. However, the amplitude of the oscillation is highly dependent on the values of 
the parameters A, B, and c. In other words, the amplitude is determined by the relation 
between the timescale of the phase transition (or the period of the oscillation) and that of 
the dust growth timescale. In the following, we discuss the relation between the the period 
of the phase transition (denoted by r tr ) and the dust growth timescale (r grow ). 

In order for the dust to grow, the dust must have enough time for the growth, i.e., 
Tacc < r tr- By using equation (H), the condition is equivalent to 

^ < r tr , (18) 

^-cold 

which reads the following inequality by using /3 grow : 

T- > T~^- (19) 

r SF Pgrow-^-cold 

Putting figrow = 100 and -X" cold = 0.5, we obtain r tr /rsF > 1/50 for the enough dust growth, 
which realizes the large amplitude of the time variation of the dust-to-gas ratio. This means 
that the timescale of the phase transition should be larger than several x 10 7 yr for the 
sufficient dust growth, since rgp ~ 3 Gyr (Kennicutt et al. 1994). Here, we note that the 
timestep in Figures 3 and 4 is c _1 in the physical unit and that cr 1 = 3 x 10 6 -3 x 10 7 yr. In 
the longer unit (c _1 = 3 x 10 7 yr; corresponding to Figures 3a and 4a), the amplitude of the 
dust-to-gas-ratio oscillation is large. On the contrary, in the shorter unit (c" 1 = 3 x 10 6 yr; 
corresponding to Figures 3b and 4b), the amplitude is much smaller because the timestep is 
shorter and thus the condition (|T9|) is difficult to be satisfied. Moreover, comparing Figures 
3a and 4a, we see that the amplitude is larger if the period is longer. 

In Figure 4a, the maximum of the dust-to-gas ratio is an order of magnitude larger than 
the minimum. Thus, the dust-to-gas ratio in spiral galaxies can show an order-of-magnitude 
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oscillation. We should note that the parameter range is reasonable for the spiral galaxies 
(Ikeuchi 1988). Thus, we have confirmed the interpretation of the scatter of the dust-to-gas 
ratio in Hirashita (1999b), who interpreted the scatter of the dust-to-gas ratio of the nearby 
spiral sample as the short-term (~ 10 7 -10 8 -yr) variation of the dust-to-gas ratio with the 
maximum-to-minimum ratio of more than 4. 

The important point is that the dust-to-gas-ratio oscillation with large amplitude 
occurs in a reasonable parameter range for spiral galaxies (Model 4a). However, the 
amplitudes are small for the other models after the stable limit cycle is attained. Thus, 
some spiral galaxies may not experience the prominent oscillation of the dust-to-gas ratio. 
But it is possible for the spiral galaxies to experience the oscillation of the dust-to-gas ratio, 
which is observed in the scatter of the dust-to-gas ratio of the spiral sample. 



4.2. Effect of Chemical Evolution 

In the discussions above, the metallicity was fixed because the timescale of the phase 
transition is much shorter than that of the chemical enrichment. However, the time 
range shown in the figures in this paper (250 timesteps) corresponds to 7.5 Gyr when the 
ctsf — 100 is adopted. Thus, we need to examine the effect of the chemical evolution here. 

To include the effect of chemical evolution, we solve equation (||) to examine the time 
variation of metallicity. If tsf is fixed for simplicity, the solution of equation (^) becomes 

X = — . (20) 

CT SF 

In other words, we simply assume a constant rate of metal enrichment to examine a effect of 
the chemical evolution qualitatively. We solve basic equations in §2.3 for the time-variable 
Xo described by equation (|20|) . As a representative case, we adopt the parameter set 
identical to Model 4a in Tablel, since the oscillation behavior of the dust-to-gas ratio is the 
most prominent (Figure 4a). The result is shown in Figure 5. We see that the qualitative 
behavior of the dust-to-gas ratio is the same as the results for constant metallicity, except 
for the gradual increase of dust-to-gas ratio. This increase is due to the metal production, 
from which dust is formed. 

Thus, two behaviors of the temporal variation is coupled: One is the gradual increase 
of the dust-to-gas ratio on the timescale of chemical enrichment, and the other is the 
short-term variation of the dust-to-gas ratio owing to the ISM phase changes. 

We note that the effect of chemical evolution is of significant importance if tsf(~ 
timescale of the chemical enrichment) is short. Especially for starburst galaxies, tsf ~ 10 8 
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yr, which is comparable to r tr (Moorwood 1996). In this case, the metal enrichment also 
proceeds from the starburst in rgp, and thus the assumption of the constant metallicity 
is never satisfied. The same discussion may also be applied to the initial burst in the 
formation epoch of galaxies, when tsf is nearly the value of starburst galaxies (~ 10 8 yr; 
the dynamical time). 

4.3. Observational Implications 

We have shown that the dust-to-gas ratio varies on a timescale of the ISM phase 
transition (~ 10 7 -10 8 yr) in spiral galaxies. The amplitude of the variation can be an order 
of magnitude (Figs. 3a-d). This confirms the previous suggestion in Hirashita (1999b) 
that the dust-to-gas ratio scatter of the nearby spiral galaxies are explained by the time 
variation. The result also suggests that the variation is partially responsible for the scatter 
of the FIR-to-optical flux ratio of spiral sample in Tomita, Tomita, & Saito (1996), since 
the dust content as well as the stellar heating (dust temperature) is responsible for the 
far-infrared luminosity (e.g., Whittet 1992, p. 170). 

We comment on the dust content in dwarf galaxies. Lisenfeld & Ferrara (1998) modeled 
the time evolution of the dust amount in star-forming dwarf galaxies (dwarf irregular 
galaxies and blue compact dwarf galaxies). Based on their model and considering the dust 
formation in cold clouds (Dwek 1998), Hirashita (1999c) presented that the dust growth 
in dwarf galaxies is negligible compared with the dust condensation from metal ejected 
by stars. This is because the metallicity of dwarf galaxies are much smaller than that of 
spiral galaxies and the collision between grains and metal atoms are not frequent enough 
for efficienct growth of the grains. Thus, the oscillation of dust-to-gas ratio through the 
variation of the dust growth efficiency is difficult in dwarf galaxies. In this case the scatter 
of the dust-to-gas ratio of the dwarf sample may be interpreted to reflect the various 
efficiency of the gas outflow (Lisenfeld & Ferrara 1998). The wind may easily blows out of 
the dwarf galaxies because of their shallow gravitational potentials (Larson 1974; De Young 
& Heckman 1994), though the distribution of dark matter largely affects the process of the 
outflow (Mac Low & Ferrara 1999; Ferrara & Tolstoy 1999). 

We should note that the timescale of the phase transition in dwarf galaxies is expected 
to be longer because of their longer cooling timescale. Since the metal-line cooling is the 
dominant cooling mechanism for the hot gas (Raymond, Cox, & Smith 1976), the small 
metallicity of dwarf galaxies, typically ~ 1/10 of the spiral galaxies, makes the cooling 
time an order of magnitude longer (> 10 8 yr). Thus, once the hot component becomes 
dominant in a dwarf galaxy, the dwarf galaxy must wait for the hot component to cool 
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down for more than 10 8 yr, during which the growth of dust grains is prevented. This as 
well as the mechanism in the previous paragraph makes the dust growth difficult. Indeed, 
the dust-to-gas ratios in dwarf galaxies are known to be much smaller than those in spiral 
galaxies (e.g., Gondhalekar et al. 1986). 

5. SUMMARY 

We have discussed the time evolution of dust-to-gas mass ratio in the context of 
multi-phase model of interstellar medium in spiral galaxies. The phase transition of 
interstellar gas occurs on a timescale of ~ 10 7 -10 8 yr in spiral galaxies (IT83). The dust 
growth rate also varies on that timescale owing to the time variation of the fractional mass 
of the cold gas through the phase transition (Hirashita 1999b; §2.1 of this paper). In order 
to examine the response of the dust-to-gas ratio to the phase transition, we modeled the 
time variation of the dust-to-gas ratio with a model of the phase transition. We adopt 
the three phase model for the interstellar gas and the Ikeuchi-Tomita model for the phase 
transition. According to the model, three types of solutions are possible: [1] all the gas is 
transformed to a hot gas; [2] a stable stationary state of three phases is realized; [3] the 
filling factors of each phase cyclically changes. For each of the three types of solutions, the 
dust-to-gas ratio behaves as follows (§3): [1] the dust is fully destroyed (the dust-to-gas ratio 
becomes 0); [2] the dust-to-gas ratio converges to a stationary state; [3] the dust-to-gas ratio 
varies cyclically in response to the phase transition. When the case [3] is applicable, the 
amplitude of the variation of the dust-to-gas ratio is large (nearly an order of magnitude) if 
the dust growth timescale is shorter than the phase transition timescale. This condition is 
satisfied within the reasonable parameter range of spiral galaxies (§4.1). 

However, the model needs to be modified for dwarf galaxies, since their small metallicity 
makes the dust growth rate in clouds small (Hirashita 1999c), and their timescale of the 
phase transition is longer owing to their longer cooling timescale. For starburst galaxies, 
the effect of chemical enrichment is important because of their short metal-enrichment 
timescale. 
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FIGURE CAPTIONS 

FIG. 1 — Upper panel: the time evolution of the the cold-gas mass filling factor (X cold ) 
for Model I. The initial condition is (X cold , X hot ) = (0.0, 0.1). The adopted values of 
the parameters are A = 1.2, B = 3.0, and ct S f = 100. The timestep is normalized by 
c' 1 ~3x 10 7 yr. Lower panel: the time evolution of the dust-to-gas ratio (V). 

FIG. 2 - - The same as Fig. 1, but for Model 2. The initial condition is 
(X co id, Xhot) — (0.0, 0.1). The adopted values of the parameters are A = 0.4, 
B = 3.0, and ctsf = 100. The timestep is normalized by c~ 1 ~ 3 x 10 7 yr. 

FIG. 3a - - The same as Fig. 1, but for Model 3a. The initial condition is 
(X co id, ^hot) — (0.0, 0.1). The adopted values of the parameters are A = 0.3, 
B = 3.0, and cr SF = 100. The timestep is normalized by c _1 ~ 3 x 10 7 yr. 

FIG. 3b - - The same as Fig. 1, but for Model 3b. The initial condition is 
(^coid, X h ot) = (0.0, 0.1). The adopted values of the parameters are A = 0.3, 
B = 3.0, and cr SF = 1000. The timestep is normalized by c _1 ~ 3 x 10 6 yr. 

FIG. 4a - - The same as Fig. 1, but for Model 4a. The initial condition is 
(X co id, Xhot) = (0.0, 0.7). The adopted values of the parameters are A = 0.3, 
B = 0.5, and cr SF = 100. The timestep is normalized by c" 1 ~ 3 x 10 7 yr. 

FIG. 4b - - The same as Fig. 1, but for Model 4b. The initial condition is 
(X co id, Xhot) = (0.0, 0.7). The adopted values of the parameters are A = 0.3, 
B = 0.5, and ctsf = 1000. The timestep is normalized by c _1 ~ 3 x 10 6 yr. 

FIG. 5 — Time evolution of the dust-to-gas ratio. The model parameters are same as 
Fig. 4a, but the effect of the chemical evolution is included. A constant rate of metal 
enrichment is assumed. 
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